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Abstract. XMM-Newton observations of the O supergiant f Orionis (09.7 lb) extend 
knowledge of its high-resolution spectrum beyond the Cvi line at 33. 7A and suggest a 
new framework for the interpretation of the X-ray spectra of single hot stars. All the lines 
are broad and asymmetric with similar velocity profiles. X-rays probably originate in the 
wind's terminal velocity regime in collisionless shocks controlled by magnetic fields rather 
than in cooling shocks in the acceleration zone. During post-shock relaxation, exchange of 
energy between ions and electrons is so slow that electron heating does not take place before 
hot gas is quenched by the majority cool gas. The observed plasma is not in equilibrium and 
the electron bremsstrahlung continuum is weak. Charge exchange, ionization and excitation 
are likely to be produced by protons. Fully thermalized post-shock velocities ensure high 
cross-sections and account for the observed line widths, with some allowance probably nec- 
essary for non-thermal particle acceleration. In general, the form of X-ray spectra in both 
single and binary stars is likely to be determined principally by the amount of post-shock 
electron heating: magnetically confined X-ray plasma in binary systems can evolve further 
towards the higher electron temperatures of equilibrium while in single stars this does not 
take place. The long mean-free path for Coulomb energy exchange between fast-moving 
ions may also inhibit the development of line-driven instabilities. 
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outflows - Shock waves 
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1. Introduction 
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The 09.7 lb supergiant £ Orionis, or Alnitak, is the easternmost star of Orion's belt and lies just 
above the Horsehead Nebula. In addition to its glorious appearance in the night sky, it is a promi- 
nent star in other respects. It is optically the brightest O star and, at X-ray fre quencies, it is the 



seven t h most luminous in absolute t erms and the fourth in apparent brightness (IBerghofer et al 



1996) 



Waldron & Cassinelli 



(2001) observed it early with the Chandra high-energy gratings to 
yield one of the first published high-resolution X-ray spectra and argued for the origin of at least 
some of the X-rays very close to the stellar surface at the base of the powerful wind that is a 
ubiquitous feature of such hot stars. The intensity ratios of the He-like triplets, which are suscep- 
tible to modification by the strong photospheric UV radiation field, were used to locate the site 
of production. 

Its X-ray spectrum has proved to be typical of those seen from O stars with very prominent 
broad lines from a variety of H-like and He-li ke io ns of common light elements and L-shell lines 
from Fe xvn in particular. 



Miller et al 



Miller (2002) commented on the unexpected 



(120021) and 

and little understood line widths observed from star to star and made general efforts to reconcile 
the data with the popular view that shocks developing from instabilities in the wind line-driving 



mecha nism are responsible for generating the X-rays. Whatever their origin, 



Feldmeier et al. 



d 19971) suggested that X-ray temperature shocks cool very rapidly via radiation losses in the 
inner wind, probably accounting for the range of temperatures of one to a few million degrees 
observed in even comparatively low-resolution ROSAT PSPC spectra. It is clear from the high- 
resolution spectra that a range of apparent temperatures is implied by the simultaneous presence 
of prominent lines such as Ovu typical of 0.8MK gas; Ovm of 2.4MK; Fexvn of 4MK; and 
Sixmof 10MK. 

XMM-Newton is ideal for observations of such hot stars. The bandwidth of the Reflection 
Grating Spectrometer (RGS) matches perfectly that part of the X-ray spectrum in which the lines 
occur, extending the Chandra High-Energy Transmission Grating (HETG) spectra beyond 25 A 
to the N vi lines near 29 A and C vi near 34A into the temperature regime below a million degrees 
where the radiative emissivity increases rapidly. The RGS is able to resolve the lines and its high 
sensitivity allows the accumulation of enough statistics to study line profiles in detail. 

The supergiant in which we are interested is the dominant component of a multiple system 



( Hu mmel et al 



2000) and is more properly known as £ Orionis Aa, though we refer to it simply 



as £ Orionis throughout. Its fainter companions are ( Orionis Ab at 0.045"and f Orionis B at 
2.4", neither of which is likely to contribute to the X-ray emission; and f Orionis C at 58", which 
would be easily resolvable by all theXMM instruments. ( Orionis B was resolved in the Chandra 
observation at an intensity of a few percent of f Orionis. Some relevant parameters which we 
have adopted for £ Orionis are listed in Table Q] Because the terminal velocity is an important 



* Based on observations obtained with XMM-Newton, an ESA science mission with instruments and 
contributions directly funded by ESA Member States and NASA 
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parameter, it is worth bearing in mind that earlier work bv lPrinia et al.1 ( 1 1990ft preferred a value 
of 1860 km s _1 , the difference probably reflecting the uncertainty in such measurements. 



Table 1. Relevant data for HD37742 09.71b f Orionis from lLamers etal.1 (119991) unless stated 
otherwise. 



visual magnitude 


V 


2.03 




effective temperature 


T e g 


30900 K 




radius 


R 


31 R 




terminal velocity 


Voo 


2100 km s" 1 




mass-loss rate 


M 


1.4 x 10~ 6 M Q yr 1 




distance 


d 


473 pc 


de Zeeuw et al. (1999) 


interstellar column density 




2.5 x 10 20 cm" 2 


Diolas & Savaee (1994) 



2. Observations 



200 lb has been in oper- 



The European Space Agency's XMM-Newton Observatory (Hansen et al. 
ation since 1999. Some of its operational and instrumental characteristics and a range of initial as- 
trophysical results were described in the special 200 1-01-01 issue of Astron omy & Astrophysics, 
including the high-resolution RGS spectrome ters by ; 



ing moderate-resolution spectrometers by 



den Herder et al 



Struderetal 



(2001) and 



(2001) and th e EPIC imag 



Turner et al 



(2001). XMM 



was used to observe ^"Ononis for nearly 12 hours in September 2002 as part of guaranteed 
time. All the i nstrum ents operate simultaneously except, in this case, the Optical Monitor (OM) 



dMason et al 



2001ft . which was not used because of the extreme optical brightness of the target. 
Some observational details are shown in Table [2] The data were reduced with standard proce- 
dures using the XMM Science Analysis System SASv6.5.0 with the calibration data available up 
to July 2005. 

The combined RGS spectrum in Fig.Q]shows that the type of strong lines seen in the HETG 
spectra below 25 A continue to the longest RGS wavelengths. The He-like O vn triplet was the 
strongest line followed by Ovm Lya; the FeXVII lines between 15 and 17 A; and the strong 
detection of C vi Lyo- near 34A. Where they overlap between about 6 and 25 A, the RGS and 
HETG spectra are essentially identical as discussed in more detail below. 

The EPIC spectra are shown in slightly different form in Fig. [2] without any effective area cor- 
rections as this becomes misleading when the energy resolution is relatively poor. Nonetheless, 
the EPIC data complement the RGS because of their higher sensitivities, as emphasized by the 
detected count rates in Table [2] and wider bandwidths. Both pn and MOS show that almost all 
( Orionis's X-rays fall in the RGS waveband. At short wavelengths beyond the strong MOS 
detections of the He-like triplets Mgxi and Sixm, apart from weak detections of Sixiv and 
S xv, there is no evidence of a hard X-ray tail. At the longest wavelengths, the counts are over- 
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Fig. 1. The RGS spectrum of ( Orionis on 2002-09-15. The spectrum is a combination of all the 
available RGS1 and RGS2 1st and 2nd order data corrected for instrumental sensitive area using 
the SAS task rgsfluxer. 

whelmingly due to redistribution from higher energies as the CCD energy resolutions deteriorate 
due to surface effects. Despite the limited ability to distinguish individual lines, the high pn 
count rate provides the most effective m eans of assessing sour c e vari ability of which there was 



none detectable during this observation. 



Berghofer & Schmitti d 19941) made a variability study 



of £ Orionis with ROSAT and found effectively no evidence for either short-term or long-term 
variability although there was one isolated episode during which the X-ray count rate increased 
by 15%. 



Table 2. XMM-Newton £ Orionis observation log showing the background-corrected count rates. 
The RGS values are the 1st order rates. 



Instrument 


Filter 


Mode Start and End Dates 


Duration(s) 


Count Rate(s~') 


RGS1 




Spec+Q 2002-09-15T13:ll:17 


2002-09-16T00:52:38 


41979 


0.346±0.003 


RGS 2 




Spec+Q 2002-09-15T13:ll:17 


2002-09- 16T00:52:38 


41979 


0.322±0.003 


MOS1 


Thick 


Small Window 2002-09- 15T 13: 12:32 


2002-09-16T00:51:00 


41728 


1.735±0.007 


MOS2 


Thick 


Timing 2002-09-15T13:12:28 


2002-09-16T00:46:41 


41473 


1.692±0.006 


pn 


Thick 


Full Frame 2002-09-15T14:04:08 


2002-09- 16T00:5 1:20 


38362 


6.659±0.014 
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Fig. 2. Spectra of ^Orionis on 2002-09-15 from the three EPIC instruments. The pn, shown 
in black, had the highest count rate mainly because it does not share its telescope module with 
an RGS. The MOS instruments have better energy resolution and were able to separate lines 
reasonably well. The agreement was good between MOS 1 and MOS2 which were operated in 
different modes. MOS1, shown in red, was in small-window imaging mode while MOS2 was in 
timing mode. 



3. Shape and strength of X-ray lines in £ Orionis 

Comparison, for example, of the Lya lines of C vi /133.734 and Nex /112.132 separated by nearly 
a factor of 3 in wavelength shows they had shapes in velocity space that were very similar as 
shown in Fig. [3] This similarity appears to apply to all the lines and we have been able to synthe- 
size a good fit to the entire observed high-resolution spectrum with the same velocity profile for 



every 



line. In order to mode l the ty pe of line asymmetries in £ Puppis discovered by 



Kahn etal 



(2001) and 



Cassinelli et al 



(2001 ). we added to XSPEcQ a local triangular line-profile model 
named TriLine characterized, as illustrated in Fig. [4] by three shape parameters: independent 
red and blue velocities where the profile goes to zero and a velocity shift of the central peak from 
the laboratory wavelength. 

The model included 85 li nes from 18 ions w ith TriLine pr ofiles, with laborat ory wave- 



lengths taken from the APED dSmith et all 1200 II) and CHIANTI (lYoung et al 



2003) line lists, 



and a thermal continuum subject to interstellar absorption fixed at the value given in Table Q] 
The 18 model ions were the H-like ions of C vi, N vn, O vm, Ne x, Mg xn and Si xiv; the He-like 



http://xspec.gsfc.nasa.gov/docs/xanadu/xspec/index.html 
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ions of N vi, O vn, Ne ix, Mg xi and Si xm; Fe xvii-xxi; and S xii-xm. The Lyman series lines were 
explicitly treated as doublets with the flux ratio fixed at 2:1. The continuum was intended to ac- 
count empirically for the combination of true electron continuum and other lines. The nature of 
the continuum is discussed further below in Section [4] 

Fits were not made to fluxed spectra like that plotted in Fig. [T]but to the detected counts 
spectra using the appropriate response matrices. Many spectral bins contain few or even zero 
counts. For such photon-limited data to which Poissonian rather than Gaussian statistics apply, 
a likelihood measure like that provided by XSPEC's C-statistic is the correct choice rather than a 
^-statistic. Errors reported in this paper were calculated for AC = 1. 

The best-fit parameters reported in Tables [3] and [4] were calculated for a combined fit to 
the spectra of RGS1 and RGS2 1st and 2nd order and the Chandra MEG and HEG ±1 spectra 
retrieved from the public archive. 

These models were a significant technical challenge with the large number of different lines 
and spectra involved, to accommodate which an expanded version of XSPEC vl 1 was built. There 
were 3472 parameters in the model which reduced to 84 for the single line profile once all the 
constraints were applied through XSPEC's Tel script capabilities. 

The RGS and HETG spectra, which were taken over two years apart, agree well with each 
other both in the shapes of the lines and in overall luminosity to within 10 or 20% as shown in 
Table|4] The single line-profile fits showed very similar red and blue velocities with a significant 
asymmetry due to a blue-shift of the line peak. The velocity errors are small by virtue of combin- 
ing several hundred independent estimates of a single profile from the lines of different spectra. 
The statistical errors take no account of any systematic differences between lines although the 
model-independent and statistical analyses reported below and the comparisons between data 
and model shown in Fig. [5] all suggest that the common profile gives a good representation. It 
establishes that the peak of the X-ray lines in £ Orionis is significantly blue-shifted by about 
-300km s _1 , similar though not as much as the shift in ( Puppis. The values of blue V and redV 
are 80% of ( Orio nis's terminal velocity of Voo = 2100km s _1 or 90% of the preferred value of 



Prinja et al.l (119901) . They are accurate enough to compare with any of the velocity parameters 
used to describe the complex shape of UV P Cygni profiles that trace cool pre-shock material. 

In order to help comparison with earlier work and allow an assessment of the significance 
of the line asymmetry, we also implemented a SkewLine model of a Gaussian line with differ- 
ent blue and red velocity half widths. Table [5] shows the results of a number of TriLine and 
SkewLine fits subject to different shift and symmetry constraints. The fits are of the same data 
and can thus be compared quantitatively with the likelihood C-statistic. The best of the eight 
fits in the Table is the blue-shifted SkewLine model which gives very similar parameters for the 
mean velocity distribution of line-emitting material as the best TriLine model. It is clear, in 
particular from comparison with the best-fit symmetrical unshifted models, that the mean line 
is very significantly skew ed and blue- shifted. The symmetrical Gaussian results agree well with 



those reported earlier by 



Miller (2002) for the Chandra data. 
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-4000 -2000 2000 4000 
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Fig. 3. The RGS C vi and MEG Ne x Lya velocity profiles of ( Orionis. The resolutions of the 
two instruments are roughly the same in velocity units at these wavelengths, allowing direct 
comparison of the two profiles. The MEG feature at about +3000 km s _1 is due to Fexvn and 
Fe xxi . 




blueV centrolV redV 

Velocity(km/s) 

Fig. 4. The TriLine model incorporated into XSPEC to describe the X-ray lines of f Orionis. A 
similar SkewLine model with unequal red and blue Gaussian widths was also used. 

Although the asymmetrical SkewLine gives a better fit, we prefer the more convenient 
TriLine parameterisation that is easier to relate to other observable quantities, velocity measure- 
ments and the line models calculated in the literature. The superiority of the SkewLine model 
may be related to ability of its wings to help account for weaker neighbouring lines not included 
in the model line list. In any case, longer-exposure data would allow a more reliable statistical 
comparison. With the single shifted, asymmetric TriLine model of Table [3] for ( Orionis, the 
line fluxes listed in Table [7] were determined, where no corrections for absorption were made. 
Here too, the measurements in common agree well with Miller's. 
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Table 3. £ Orionis best-fit emission-line velocity parameters and underlying continuum for the 
simultaneous multiple common-profile TriLine XSPEC model fit to the XMM-Newton RGS 1st 
and 2nd order spectra of 2002-09-15 and Chandra MEG and HEG ±1 order spectra of 2000-04- 
08. 

TriLine velocities 
blueV -1642±22 kms" 1 

centralV -302±29 km s" 1 

redV +1646+26 km s" 1 

bremsstrahlung continuum 
N H 2.5 x 10 20 cm" 2 

kT 0.494+0.007 keV 

normalisation 5.66+0.14 x 10" 3 
C-statistic=22862.4 using 26281 PHA bins 



Table 4. Best-fit normalisation parameters for the model of Table[3] TriLine model fit to high- 
resolution spectra of £ Orionis taken with XMM-Newton on 2002-09-15 and Chandra on 2000- 
04-08. 



Instrument 


Order 


Counts 


Relative Normalisation 


RGS1 


-1 


15136 


1.02+0.01 


RGS 2 


-1 


14133 


1 


RGS1 


-2 


2771 


1.13+0.01 


RGS2 


-2 


2507 


1.00+0.01 


MEG 


-1 


5556 


1.11+0.02 


MEG 


+ 1 


3568 


1.03+0.02 


HEG 


-1 


1129 


1.18+0.03 


HEG 


+ 1 


839 


1.12+0.03 
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Table 5. Best-fit emission-line velocity parameters in km s _1 for variously constrained TriLine 
and blue and red half-Gaussian SkewLine model fits to the grating spectra of £ Orionis from the 
XMM-Newton RGS and the Chandra MEG and HEG. 

blueV centralV redV C-stat constraint 

1642±22 -302+29 +1646±26 22864.4 Best fit of TablefJ] 

1641+15 -303±26 +1641+15 22865.9 |redV|=|blueV| 

1774+21 +1492+21 22973.1 centralV=0 

1684+15 +1684+15 23053.2 centralV=0,|redV|=|blueV| 

blueW centralV redW C-stat constraint 

SkewLine 650+25 -318+29 991+27 22816.7 

SkewLine 811+9 -126+10 811+9 22879.4 redW=blueW 

SkewLine 880+12 738+12 22947.6 centralV=0 

SkewLine 841+10 841+10 23038.4 centralV=0, redW=blueW 



TriLine 
TriLine 
TriLine 
TriLine 



10 



Pollock: A new O-star X-ray paradigm 




Fig. 5. Selected comparisons of high-resolution X-ray data with the best-fit single mean TriLine 
profile model used for all the lines with the 1st order RGS1 in red, RGS2 in blue and the com- 
bined MEG 1st order in green. The common model is shown by the smooth black line. The O vm 
Lya line was affected by bad pixels in RGS 1 . 



Pollock: A new O-star X-ray paradigm 1 1 

3.1. The composite model-independent line profile 

The important model-independent question arises of whether any X-ray gas was observed at ve- 
locities greater than Voo. We have calculated the distribution of all the detected events in velocity 
space by stacking all the individual line velocity distributions. Any continuum or background 
photons should form a smooth background in the composite velocity space. The resultant pro- 
files from both XMM and Chandra data are shown in Fig. [6] It should be emphasized that this 
plot has limitations in that it includes the instrumental resolutions and neglects sensitivity vari- 
ations, such as bad pixels, that are properly taken into account in the response matrices; the 
RGS 1st order profile, for example, was affected by CCD columns missing near the centre of the 
O vin Lya /118.967. As these defects have nothing to do with the distribution of moving gas in 
f Orionis, we should be able assess reasonably objectively from the ensemble of lines the range 
of velocities present. Taking into account the differing resolving power in the four spectra, they 
agree well with each other and with the TriLine velocities reported in Table [3] The blue wing 
is sharp and joins the weak underlying continuum before the terminal velocity is reached. The 
red wing, on the other hand, flattens above about +1000 km s _1 because of the strength of the 
He-like triplet intercombination lines. It is quite safe to conclude from Fig. [6] that there was no 
X-ray emitting gas moving at or above the wind terminal velocity. 



3.2. An assessment of line-profile variations 



The sharpness of the composite line profile is consistent with the idea suggested by Fig.[3]that all 
ions have much the same shape. In order to check this more rigorously, we have estimated indi- 
vidual best-fit line profiles for the most prominent ions as reported in Table [6] As well as show- 
ing the individual ion TriLine parameters, the statistical significance of the differences may be 
judged using AC, the decrease in the value of the C-statistic gained by freeing an ion's TriLine 
parameters. Under the null hypothesis of a common line profile, AC should be distributed as x\ 
for the 3 extra degrees of freedom. The variations between ions are small, amounting to less than 
about 200 km s _1 from one ion to another, but probably real. Nevertheless, the assumption of a 
single profile looks like a sensible working h y pothe sis. 



It has been argued by both 



Kahn etal 



(2001) and 



Cassinelli et al 



(2001) that the Nvn 



Lya line in f Puppis has a much wider flat-topped profile, significantly different from the other 
lines, so that a single shape does not apply. However, the RGS data of ( Puppis suggest that the 
N vn line might not be unusually wide if account is taken of blending with N vi He/3: the N vi 
He-like triplet is strong and in common with N vn Ly/3 and C vi Lya is not obviously unusually 
broad. We intend to discuss £ Puppis in detail in a future article. In £ Orionis, the evidence from 
Tableland FigQis that the longest-wavelength C vi line is, if anything, narrower than average 
rather than broader. 
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Table 6. Best-fit TriLine model velocity parameters for individual prominent ions in the X-ray 
spectrum of f Orionis compared with the single common line model. The quoted lcr errors were 
estimated using the XSPEO error 1 . command 



ion 


blueV 


centralV 


redV 


C-stat 


AC 




(kms- 1 ) 


(km s" 1 ) 


(km s- 1 ) 






All ions 


-1642+22 


-302+29 


+1646+26 


22862.4 




Cvi 


-1719+135 


-270+137 


+ 1403+130 


22855.6 


6.8 


Nvn 


-1705±55 


-111+240 


+1065+324 


22859.8 


2.6 


Ovn 


-1349±72 


-509+83 


+ 1835+76 


22844.9 


17.5 


O vin 


-1656±43 


-258+59 


+ 1482+46 


22838.7 


23.7 


Fe xvn 


-1647±42 


-421+62 


+ 1845+59 


22848.9 


13.5 


Ne ix 


-1653±77 


-211+79 


+1508+102 


22859.9 


2.5 


Nex 


-1730±106 


-327+126 


+ 1735+103 


22861.0 


1.4 


Mg xi 


-1536±153 


+33+121 


+ 1487+144 


22850.1 


12.3 


Mg xn 


-721±365 


-662+227 


+ 1514+471 


22855.1 


7.3 


Si xiii 


-1449±124 


+48+240 


+ 1421 + 188 


22856.4 


6.0 
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Fig. 6. Composite X-ray line profiles of £ Orionis calculated from XMM-Newton RGS and 
Chandra HETG grating data. The instrumental resolution improves down the page so that the 
observed profile becomes a better approximation to the actual velocity dispersion of X-ray emit- 
ting material. The profile at v > +1000 km s _1 is distorted by the intercombination lines of the 
He-like triplets, f Orionis's terminal velocity is Voo = 2100km s _1 . 
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Fig. 7. The C vi Lya line of £ Orionis observed with the XMM-Newton RGS shown by the his- 
togram compared with the mean TriLine model profile of the joint fit of Table[3]to all the lines 
in the spectrum. 
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4. The X-ray continuum in f Orionis 

Although a semi-empirical continuum component was added to the line models described above, 
the real nature and strength of the continuum is an important matter as it is a direct diagnostic of 
the density and temperature of plasma electrons. In common with other O-stars, the continuum is 
weak. Not only are lines responsible for the most obvious features in the spectrum in Fig.Q] but 
increasing numbers of weaker lines seem to make up most of the emission observed. Fig.[8]shows 
details of the spectrum between 17 and 19 A. The low flux between about 17.8 and 18.5 A, that 
contrasts to the plateau immediately bluewards, is common to the handful of O stars observed so 
far but still shows features that could plausibly be identified with upper level transitions of Fe xvm 
for which no observed wavelengths are available in the line lists. The flux reached similarly low 
but again not especially smooth values between 20 and 21 A, where plausible identifications are 
with lines of N vn and Ca xvn. 

The high-resolution X-ray spectra of well known relatively high-temperature active cool 
stars, such as AB Dor and HR 1099, show obvi ous strong continuum emission, although this 



is not the case with cooler stars such as Procyon (IRaassen et al 



2002), where the continuum is 



notably weak. We considered global empirical line+continuum models for £ Orionis with more 
than the 85 lines already included but these soon became impractical as the number of weak lines 
increases exponentially and the wavelengths are uncertain. If a continuum is present in ( Orionis, 
its flux near 18A is less than about 5 x 10 _5 cm _2 s _1 A -1 . A long exposure would help to constrain 
it with more confidence. 



16 
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XMM-Newton RGS £ Orionis 

0.0010 1 I' 1 III f 1 1 I 1 1 II 




Wavelength(A) 

Fig. 8. The XMM-Newton RGS spectrum of £ Orionis between the lines of 
Fexvn/U17.053, 17.098 and Ovu ^18.627 and Oviii XU8.967, 18.972. Some of the weak 
lines between about 17.8 and 18.4 A are probably transitions to the excited 2s2p 6 level of 
Fe xvm, which appear in the ATOMDB and CHIANTI line lists with uncertain wavelengths. 
The level of the continuum looks to be less than about 5 x 10~ 5 crrr 2 s~' A -1 in this part of the 
spectrum. 
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5. Origin of the X-rays from C Orionis 

The measured fluxes and velocity widths of the lines in the X-ray spectrum of £ Orionis are 
plotted in Fig. [9] The fluxes have been converted to emission measures usi ng the ATOMDeP 



maximum emissivi ties, similar to the plots against temperature for £ Puppis by 



Kahn etal 



(2001) 



and 5 Orionis by Mil ler et alJ (12002). In common with those and other stars, the X-ray emission 
measures are three or four orders of magnitude below that available in the wind as a whole and 
show that only a small fraction of the wind material is involved. 

The remarkable thing about Fig.[9]is the lack of any substantial variation with wavelength of 
either emission measure or velocity width, in contrast to the rough expectations of the instability- 
driven shock model also plotted there. These have been calculated for the parameters of f Orionis 
in Table Q] assuming a ve locity law v(r) = Voa(l — RJ rf with f3 = 0.8. In the development of 



standard model following 



Waldr on & Cassinellil (1200 ll) 's first discussion of ( Orionis's Chandra 



grating spectrum, X-rays at shorter wavelengths are thought to emerge from deeper into the 
slower, denser material of the acceleration zone of the wind as the opacity decreases, with the 
corresponding expected increase in the emission measure plotted in Yet the measurements show 
little or no evidence of this. Rather they give quantitative expression to the impression of FigQ] 
that the spectrum is more-or-less optically thin. If this is the case, then the suggestion is that 
X-rays are not from the dense central regions of the acceleration zone but from further out in 
the terminal velocity region, although a more precise location is not known. How this might be 
reconciled with the apparently firm conclusions concerning the location of the X-ray emission 
deduced from analysis of the He-like triplets is one of the topics discussed below. 

In the following sections, an attempt is made to discuss in some detail the feasibility and 
consequences of the possibility that X-rays arise in the wind's terminal velocity regime. It should 
be stressed that, in common with any ideas, success or failure should be judged by three things: 
consistency with the law of physics; ability to reproduce current data, such as those in Fig- [9j and 
ability to allow predictions to be made that can be tested against future observations. 

To begin with, it is tempting to interpret the simultaneous presence in the spectrum of lines 
of a range of ionization potentials as proof of the presence of a range of temperatures, lending 
support to the proposal that cooling rather than isothermal shocks are responsible for the X-ray 
emission. For £ Orionis, using the well-known methods offered, for example, by XSPEC or SPE3U 
it is possible to make a good synthesis of the spectrum with a multi-temperature model involving 
gas in collisional ionization equilibrium between about 1MK and 7MK involving either a limited 
number of discrete temperature components or a continuous DEM distribution, like that sug- 
gested by Fig. [9] However, before accepting at face value the tempera tures that this would imply , 



it is worth considering the general physical principles discussed by 



Zel'dovich & Raizeri (2002) 



that govern shocks and how these might affect their development in the particular circumstances 



2 http://cxc.harvard.edu/atomdb/ 

3 http://www.sron.nl/divisions/hea/spex/ 
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Fig. 9. Observed properties of the X-ray lines of ( Orionis. In the top panel are plotted the line 
emission measures, EM - 4-nd 2 (f/e mix ), calculated from the measured fluxes, /, of Table [7] 
the maximum solar-abundance APED emissivities, e max , and the distance, d, to the star. The 
small error bars reflect the high precision with which the fluxes are known. Below are plotted 
the measured redV, centralV and blueV parameters of the TriLine models of individual line 
profiles from Tabled with the corresponding mean values shown by the horizontal lines. Also 
shown as solid lines are two illustrative quantities calculated for R T= i(A), the radius of the t = 1 
X-ray absorption optical depth for a /3 = 0.8 velocity law. At the top is EM(> R T= i{A)), the 
emission measure outside R T= \(A), scaled down to overlap the data. Below is +v(R T= i(A)), the 
range of velocities offered by the accelerating wind at that radius. 
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of a stellar wind , especially if the sh ocks occur in the relatively low density terminal velocity 



region. Recently 



Pollock et al 



(2005) did this type of analysis for the binary-system colliding- 
wind shocks in WR 140. Similar considerations for single hot stars lead to some interesting 
conclusions. 
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5.1. The nature of shocks in O-star winds 

An O-star wind is a plasma flow in which particle interactions are long-range Coulomb colli- 
sions. These have a key role in the wind acceleration mechanism in normal circumstances as 
Coulomb coupling is responsible for redis tribution of momentum from the small minority of 
UV-driven ions to the rest of the flow (e.g. '. 



Kudritzki & Puis, 



2000). In more extreme conditions, 



ZeFdovich & Raizerl d2002l) described the particular properties of plasma shock waves that result 
from the slow character of the Coulomb energy exchange among ions and even slower between 
ions and electrons. Because most of the kinetic energy of the flow is in the ions, the most im- 
portant quantity that determines shock development is the ion-ion collisional mean-free path, 
which can be estimated as the distance required to deflect by 90" through successive Coulomb 
collisions a proton or H e nu cleus moving with the velocity, v = vs x 1000km s _1 . Following 



Draine & M cKee ( 1993) after 



Spitzerl d 19621) . the mean-free-path is 



h-i ~ 7.0 x 10 18 v 8 4 M cm 



(1) 



where n\ is the ion density. Fig. [TU] shows how the ion-ion collisional mean-free path varies 
throughout the wind of £ Orionis. 

Close to the photosphere, where the density is high and the velocity low, is small, but 
then increases rapidly with the increasing velocity and falling density. At about 3R* above the 
stellar surface in the acceleration zone, « Q.lR t while at 10R t , x R t . The implications 
are profound. If strong shock discontinuities are to develop at all in an O-star wind then some 
other dissipation mechani sm than collisions must c ome into operation: they must be collisionless 



shocks. As discussed by iDraine & McKed (11993b . for example, collisionless shocks probably 



account for the majority of shocks observed in cosmic plasmas, which are of predominantly 
low density and thus of long collisional mean-free-path. Familiar examples are the interaction 
between the solar wind and the earth's magnetosphere and the blast wave of a SNR sweeping 
through the interstel lar medium; more relevant might be the colliding-wind shocks of WR 140 



dPollock et al 



20051) . In all these cases, shock dissipation occurs though collective plasma pro- 



cesses involving a magnetic field in which a small characteristic dissipation length is set by the 
ion Larmor radius. The development of shocks in single stars and the production of X-rays is 
thus likely to be a phenomenon controlled by magnetic fields at large within the wind. 
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Fig. 10. The mean free path for ion-ion Coulomb collisions in the wind of £ Orionis. 
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5.2. Post-shock relaxation and the electron temperature 



As described bv lZei'dovich & Raizerl (120021) . a property of plasma shock transitions is that ions 
are heated to high temperatures while electrons initially remain relatively cold. This is certain 
for shocks in media dense enough for Coulomb collisional dissipation to op erate but is proba - 



Rakowski (2006) 



bly also largely true in collisionless shocks, according to the recent review by '. 
Although the amount of electron shock heating by plasma turbulence above the minimum pro- 
vided by adiabatic heating remains a matter of speculation, Rakowski showed that the evidence 
from SNR shocks is that electron heating decreases with shock velocity and should be small for 
the few thousand km s typical of stellar- wind terminal velocities. It is in the subsequent flow 
through the post-shock relaxation layer that energy is transferred from ions to electrons as the 
plasma moves towards equilibrium. Naturally enough, for the use of equilibrium plasma models 
to be justified in the post-shock gas, equilibrium between ions and electrons needs to have been 
established. 

Assuming post-shock bulk energy redistribution takes place through Coulomb collisions, two 
stages are expected, first amongst ions and then between ions and electrons. During both these 
stages, the ionization state will also be moving from that of the pre-shock gas towards its relaxed 
state. Post-shock equilibrium would apply once all these various stages have run to completion. 
We suggest that the time-scales are too l ong fo r this to happen. 



Spitzerl dl962l) . iKrall & Trivelpiecd d 19861 section 6.4) and 



Zel'dovich & Raizerl (120021) . for 



example, showed that equilibration between ions and electrons is much slower than between 
electrons or ions among themselves because of the large difference in mass between ions and 
electrons. I Spitzerl (1 19621) estimated about a factor of 50 slower than ion-ion equilibration from 
the square root of the mass ratio for modest departures from equipartition. 

Immediately behind a stellar wind shock, conditions are far from equilibrium. iBrown et al 



( 120051) have made rigorous calculations of the rate of temperature equilibration, incorporating 
simultaneous Coulomb scattering on short scales and collective plasma heating on large scales, 
for arbitrary ion and electron temperatures, T, and T e . Energy transfer is governed by the term 



2„2 



d(Ti - T e )/dt oc K fK e , , 

(A-m,+ fame?! 1 

where i<f = y8,e^n, and = This equation shows that in the very earliest stages of re- 

laxation, when ion and electron random velocities are similar and T- t » T e , the electrons will 
initially heat up quickly in a time proportional to (m,m e ) intermediate between those of electrons 
and ions among themselves. As the electrons heat up, the fraction of energy (m e jmi) transfered 



in a collision has less effect and the rate slows down to the factor yfmjmi times ion-ion rate 
calculated by Spitzer. 

Thus, the length scale for post-shock electron heating in the wind of an O star compared to the 
ion-ion mean-free-path of Fig[l0]is likely to be many stellar radii in which case the hot gas may 
not be able to survive long enough for equilibrium to be established. The low X-ray emission 
measures of O-stars show that only a small fraction of the wind material is involved, the vast 
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bulk of which is cold: shocked gas will soon get mixed back with cool material and disappear 
from view. In this time, there will be little chance for electrons to be heated enough to be able to 
contribute to the observed X-ray emission by either ionization or bremsstrahlung continuum. A 
scarcity of shock-heated electrons in the winds of O stars would give rise to a weak continuum 
and require an alternative means of producing the high levels of ionization. 



5.3. X-ray ionization by ion-ion interactions in O-star winds 

In the supposed absence of shock-heated electrons, what is the origin of the ion species un- 
ambiguously identified in the X-ray spectrum? Ions in general and protons in particular in the 
imme diate post-shock gas are probably responsible, as anticipated in general terms by 



Mewe 



(1999, p. 148). Through the initial shock transition, the ionization balance is not changed al- 
though ions characteristic of the partially ionized cool wind suddenly find themselves in a hostile 
environment: encounters take place with other ions at relative velocities caused by randomization 
of the bulk velocity of the pre-shock gas of about 2000 km s _1 in the terminal velocity regime 
of f Orionis. Protons of such velocities are a similarly effective agent for ionization and exci- 
tation as electrons of a few keV because the cross-sections depend on the relative velocity of 
the incident ionizing particle and the velocity of the bound electron, whose order of magnitude 
is fixed by the Bohr velocity VBohr = 2188 km s _1 . We suggest that it is the coincidence of this 
microscopic atomic value of v'Bohr - widely used as the atomic unit of velocity - and the macro- 
scopic terminal velocities of ( Orionis and other O stars, that is the basic physical reason for the 
production of X-rays in hot stars. 

The relaxing material contains no neutral atoms but is still far from the fully-ionized state 
discussed by Spitzer and others, so that the same ion-ion interactions involved in energy ex- 
change via Coulomb collisions are likely simultaneously to affect the remaining bound elec- 
trons, with possible observational consequences. The ion-ion Coulomb cross-section of 1 ~ 
7 x 10" 21 cm 2 estimated from equation ([TJ, which determines the rate of energy exchange in elas- 
tic collisions between ions in the post-shock gas, serves as a benchmark in the discussion below 
against which the effectiveness may be judged of inelastic processes that also occur during re- 
laxation. High-velocity encounters between ions have maximum inelastic cross-sections of order 
10" 16 cm 2 that exceed the Coulomb cross-section by orders-of-magnitude and thus might be 
expected to dominate post-shock relaxation as discussed in more quantitative terms below. 

Ion-ion interactions that involve bound electrons are of three types, namely ionization, ex- 
citation and charge exchange. Ionization and charge exchange are both mechanisms of electron 
loss that lead to changes in the ionization state and so could be responsible for the highly-charged 
ions observed in the X-ray spectrum. Excitation and sometimes charge exchange leave ions in 
excited states from which decay leads to observable radiation including X-rays. A typical recom- 
bination cross-section is of order 10~ 21 cm 2 and too small to compete with radiative de-excitation 
in X-ray production. 
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5.3.1 . Ion-ion ionization 



Ion-atom ionization has been discussed in detail by lKaganovich et al.l (I2006I) . whose scaling law 
for the ionization cross-section by a fully-stripped projectile ion of charge Z p and velocity v of 
one of the N„i electrons bound in a quantum state (n, I) of ionization potential /„/ to a target 
nucleus of charge Zj is given by 



cr ion (v,I nl ,Z p )=na 2 



(Zp/Zj + l) 



(3) 



v„i^Z p /Z T + 1 J 

where atomic units for cross-section, velocity and energy are defined by the Bohr cross-section 
na^ = 0.880 X 10~ 16 cm 2 ; vo = VBohr = 2188 km/s; and Eq = m e vi = 27.2 eV, respectively. This 
equation shows in particular how ion-ion ionization cross-sections depend sensitively through the 
scaling function G(v/v max ) on the relative velocity of the incident ionizing particle and the bound 
electron. The maximum cross-section, <x max , is given by G(l) ~ 0.8 at v = v max = v„i **jZ p /Z T + 1. 



The velocity of a bound electron is roughly v„i/vq ~ Zr/n ~ -\J2I„i/Eo- In the immediate 
post-shock gas of a stellar wind, the relative velocity scale of ion-ion encounters is fixed by the 
randomized terminal velocity v M , though the Maxwellian distributions will en sure that a fraction 



of enc ounters also take place at relative velocities up to a few times this value. iKaganovich et al 



(2006) showed that their scaling law 



G(x)= 6XP( } IX - [1.26 + 0.283 ln(2x 2 + 25)], 



(4) 



which falls away from the maximum towards the lower velocities of interest here, reproduces well 
a variety of observational data for x > 0.5, at which point the ionization cross-section has fallen to 
20% of its maximum value. At low velocities, measurements become difficult because, to quote 
Kaganovich et al., the ionization cross-section is completely dominated by charge exchange, 
whose cross-section is comparable to cr max . 



5.3.2. Ion-ion charge exchange 



Charge exchange between ions is likely to be of fundamental importance in the immediate post- 
shock gas in the wind of a hot star and be r esponsible for production of many highly-charged ions. 



Following Bransden & McDowelll (11992. pp. 11-15), a simple theoretical estimate for the low- 
velocity charge exchange cross-section between a bare nucleus projectile of charge Z p interacting 
with a target nucleus of charge Zj with a single bound electron of principal quantum number n 



may be written 



cr CE (Z p ,Z T ) = na^ 



'2n 2 (Z T +2yfZpZ;y 



(5) 



independent of velocity for v/v„; < (Z p jAZj)i . This expression reproduces experimental data 
well for electron capture by ions from atomic hydrogen and if applied to collisions between 
protons (Z p = 1) and, say, 7+ (Zj = 8) ions, predicts a cross-section of 1.6 x 10~ 17 crrT 2 . 
Crossed-beam measurements of collisions between other ions have been reported, for example, 
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Skieraetal 



(2001), 



von Diemar et al 



<|200l|)and 



Brauning et al, 



(2005) and show similar or 



by 

higher cross-sections, including more than 10~ 16 cm~ 2 for the production of N 5+ from charge 
exchange with a-particles: this is interesting because N vi lines appear in the X-ray spectrum of 
£ Orionis. 



5.3.3. Implications for post-shock relaxation and the X-ray spectrum 

These high cross-sections, which overwhelm Coulomb collisions, guarantee that the first thing to 
happen to ions in the post-shock relaxation layer concerns their outer bound electrons with charge 
exchange and ionization leading to rapid readjustment of the ionization balance. This even takes 
place faster than relaxation among different ions to a common te mperature. This conclusio n is 



Hughes & Helf and ( 1985J) for 



probably of more general significance. In SNRs, as calculated by|] 
example, it has often been assumed that the post-shock ionization balance changes in response to 
hot electrons, whereas particul arly for the shocks velocities of a several thousand km s -1 of some 



of the remnants discussed by 



Rakowskil (120061) . ion-ion interactions are als o important, as in 



the ch arge-exchange layer immediately behind some SNR shocks discussed by 
Jl993l p. 403), for example. 



Draine & McKee 



For O-star winds, we propose that it is also ion-ion interactions and notably charge exchange 
during the initial phase of post-shock relaxation that give the observed X-rays. While there seems 
little doubt that the necessary ions can be produced, the generation of the X-ray spectrum is more 
uncertain. This may be due to excitation by ions, mainly protons; charge exchange into excited 
states in the hot gas; or charge exchange in eventual encounters with the much less tightly bound 
electrons in the plentiful supply of much cooler ions in nearby unshocked material that makes 
up the great majority of the wind. This final possibility seems worthy of future work. If proton 
ionization is significant, while it is not expected that the great majority of free electrons will be 
heated much immediately behind the shock, a small population of hot secondary electrons would 
be ejected, giving some scope for a bremsstrahlung continuum. 



5.4. The velocity profile of the X-ray lines 



Most of the efforts made so far t o account for t he sha pe o f the X-ray lines in ho t stars, such as 



those of llgnace & Gayleyl (2002) 



Kramer et al 



420031) and 



Oskinova et al 



(2004), have assumed 



that the X-ray emitting gas is moving with the majority cool gas. By doing so and appealing to 
the velocity law of the cool gas, attempts have been made to constrain the location of the hot gas. 
In particular, it has regularly been argued that the red wing of the X-ray lines arises in material 
in the wind on the far side of the star flowing away from the observer. This simple assumption 
is hard to justify, if only because much of the ordered motion of the wind needs to be converted 
to heat for X-rays to be observed at all and, for CIE models, half the energy has to be transfered 
from ions to electrons for excitation to occur. 
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For an ideal gas of mean particle mass m flowing at velocity v into a strong stationary 
shock, the jump conditions determine the equilibrium temperature and velocity of the outflow- 
ing shocked gas to be kTs/m = (3/16)v 2 and vs = v/4. For material of solar composition, 
kTs ~ 1.2v? keV, where v = vg x 1000 km so that the observation of X-ray gas in equilib- 
rium near the 0.5 keV of Table [3] would imply the complete dissipation of about 650 km s _1 of 
directed kinetic energy. 

Out of equilibrium, it is more likely that the observed line velocity profiles mostly reflect 
instead the line-of-sight component of the thermalized motion of ions in the immediate post- 
shock gas. For a Maxwellian distribution, 



HWHM(vx) = V21n2(Jtr s /m) = ( V61n2/4)v ~ 0.51 v. (6) 

The similar observed widths of different ions further implies that even equilibration between ions 
does not take place. 

In addition to this thermal component, the bulk flow away from the shoc k at vs = v/4 will 



Ignace & Gavlevl (120021) made for material 



also contribute to the line profile. The calculations that|] 
moving at constant radial velocity in the terminal velocity flow, may equally well be applied to 
this slower-moving material. The resultant, roughly triangular, line profiles show blue-shifted 
peaks at 10 or 20% of the outflow velocity, depending on the optical depth, and a half width of 
about 60%. The resultant combination of thermal and bulk motions should then shift the line to 
the blue by a fraction of the terminal velocity with a width mostly determined by the greater 
random motion. 

The observed lines in ( Orionis are roughly consistent with this scheme, showing blue-shifted 
peaks of about v x ll and half-widths of about 80% of the value of 0.51voo expected for complete 
thermalization of gas moving at the terminal velocity. While it is conceivable that the X-ray 
emission is picking out the small amount of gas in the wind-acceleration zone travelling at this 
fraction of the terminal velocity, we consider this unlikely, not least because of the apparent 
lack of X-ray absorption.. Two further possibilities come to mind: either that the shocks are not 
stationary in the wind but are moving outwards at a few hundred km s" 1 and are completely 
thermalized; or that the shocks are stationary and occur throughout the terminal velocity regime 
and the energy balance involves other channels than merely the production of heat. As non- 
thermal particle acceleration is a common feature of collisionless shocks, such as those in SNRs 
and those probably in WR 140, we consider this possibility the more promising, although free- 
free emission and absorption by the cool wind is expected to be strong enough to prevent any 
non-thermal radio radiation from being observed. 



5.5. Absence of equilibrium in the post-shock plasma and the He-like triplets 

If O-star X-rays are produced by terminal-velocity shocks in the wind, then the sluggishness of 
energy exchange throughout the post-shock gas will prevent the establishment of equilibrium 
before the hot gas is overwhelmed by cool material. Therefore, the use of coronal equilibrium 
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models to analyse O-star X-ray spectra is unlikely to yield any physically meaningful information 
about conditions in the hot gas, though it is always possible to derive equivalent equilibrium 
temperatures like those quoted above for f Orionis and elsewhere for other stars. The X-ray 
spectrum is more likely to reflect the individual cross-sections of the three important processes 
between ions discussed above of charge exchange, ionization and excitation. 

An argument in su pport of the instability-driven shock model has been its success following 



Kahn etal 



(2001J)and 



Wald ron & Cassine lli (2001) in accounting for the characteristic patterns 



of the fir lines of the He-like triplets: in contrast to most other high-resolution X-ray spectra, 
the O-star intercombination lines are particularly strong, showing the effects of photospheric UV 
pumping of the 3 Si upper level of the forbidden line to 3 Pj 2 upper levels of the intercombination 



line. In ( Puppis, 



Kahn etal 



( 200 1 ) showed that there is enough radiation that the photoexcitation 
rate exceeds the intercombination line decay rate with a radius that varies from about 3R* for 
Si xm to 200/?* for N vi, so that the UV affects an enormous volume of the wind. 



(2001) 



Calculations of the expected line ratios such a UV field, such as those o flPorquet et al. 
have been based on otherwise equilibrium conditions. In the rarefied gas immediately behind 
the collisionless shock in an O-star wind imagined above, apart from the intense UV radiation 
environment, other conditions are expected to be quite different with the He-like ratios subject to 
a variety of considerations, namely: 

- a plasma out of equilibrium 

- a scarcity of hot electrons 

- energy concentrated in interacting ions 

- widespread charge exchange 

- ionization by non-thermal particles 

Hot protons and cool electrons are expected to have similar random velocity distributions 
and thus both contribute to collisional coupling of the He-like triplet levels. Proton excitation 
leads to different line ratios because of the absence of resonances and exchange terms (e.g. 

p. 41). 



Shevelko & Vainshtein, 



1993 



Charge exchange is likely to be as important in stellar winds as it is for the generation 
of cometary X-rays by charge exchange of neutral material with highly-charged ions in the 
solar wind. In laborator y measurements designed to simulate cometry charge exchange by 



Beiersdorfer et al 



(2003), which are at very low velocities on the atomic scale, the forbidden line 
was strongly enhanced. In the crossed-beam measurements discussed above at higher relative ve- 
locities more relevant to hot stars, though still low on the atomic scale, the total cross-sections 
are dominated by charge exchange. 

Although, following the suggestion by 



Pollockl (11987), most authors have cited inverse- 



Compton emission as a non-thermal signature, ionization losses are a more important cooling 
mechanism of mildly relativistic particles. This may be an opportunity to observe, at least in- 
directly, the effects of cosmic -ray acceleration in single stellar winds through the detection of 
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satellite lines following inner-shell ionization (e.g. iMewe & Schriivei 



19781) . which may also 



cause the unusually strong H e-like forbidden lines in the colliding-wind spectra 77 Carinae and 



WR 140 (Pollock etal 



20051) . Incidentally, there may be some evidence of satellite lines in the 



O vii triplet of £ Orionis in which the intercombination line appears significantly blue-shifted 
from the wavelength expected from the ensemble of other lines by -13 ± 3 mA. This partly ac- 
counts for the extra blue shift of O vn in Table [6] 

In summary, some theoretical effort will be required to calculate the X-ray spectrum in gen- 
eral and the He-like line ratios in particular to be expected from shocked ions in an O-star wind 
before proper conclusions can be drawn about the effects of UV radiation and the implications 
for the location of the X-ray emitting material. 



6. X-ray absorption in O-star winds 



In spit e of the high optic al depths expected on the basis of supposedly well-determined mass-loss 



rates, 



Kahn etal 



(1200 ll) showed the lack of any detectable absorption edges i n 4" Puppis. Even 



Paerels & Kahn, 



curso ry inspection of £ Orionis's RGS spectrum or those of other O-stars (e.g. 
2003) shows no obvious sign of absorption, as emphasized by the strength of C vi Lya. 

An important part of attempts to reconcile the observational data with an origin deep in the 
wind's acceleration zone within a few stellar radii of the surf ace has involved seeking mea ns of 



reducing the X-ray op acity. From a theoretical point- of-view, 



Cassinelli et al 



(2001) used the reduced opacities of 



Waldron & Cassinelli (2001) and 



Waldron et al 



1998) from self-consistent 



calculations of the ionization b alance to be expe cted in the presen ce of UV and X-radiation 



(IMacFarlane et al. 



1993 



1994). Independently. 



Kramer et al 



(2003), for example, appealed to 



significantly reduced attenuation in their models of the £ Puppis X-ray line profiles. 

It might have been expected that the high UV outer-shell optical depths necessary for the 
line driving would guarantee a correspondingly high X-ray inner-shell optical depth. It has been 
known for many years that this is indeed the case, as demonstrated by the bright neutro n star hard 



Sako etal. 



X-ray sources powered by accretion from an otherwise normal hot star. Vela X-l (e.g. 
( 1999) and references cited there) is a relevant example as the primary star HD 7758 1 is a B0.5Iab 
star of similar spectral type to ( Orionis. The neutron star orbits at about 0.78/?* above the surface 
of the supergiant, in what would otherwise be the strongest region of the acceleration zone. It is 
a simple observational fact that the neutron-star secondary's X-rays are subject to heavy absorp- 
tion. While this varies through the orbit, reaching values of several times 10 23 cirT 2 , much of it is 
due to the rearrangement of wind material into an accretion wake that trails the neutron star. The 
X-ray absorption, though occasionally extremely high, is well modelled by cold material. 

Immediately after eclipse, on the other hand, when the line-of-sight passes only through 
the stellar wind, a minimum column density of about 2 x 10 22 otT 2 is observed.. This value 
is reasonably close to that expected if £ Orionis were to take the place of HD 77581. For a 
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spherically symmetric wind with a velocity law, the integrated column density between a point 
distance R from the surface of the star along the line-of-sight scales as 

N R (R, <f>, i) ~ (M/fwipV^RtXy/ siny)H(R,R m ,0) cnT 2 (7) 

where /u is the mean atomic number, m p the proton mass and v(r) = v 0O (l - R t /rf. The column 
density is the product of a column density characteristic of the star, N* = (M/fmipVccR*), and 
two scaling factors: a geometrical factor (yj siny), where cosy = cos sin 6 and <f>,6 are the 
azimuth and inclination angles between the radial direction through R and the line-of-sight; and 
an analytical density -related factor H(R,R t ,fi), which exceeds unity for points in the densest 
parts of the wind and equals (R t /R) far from the stellar surface. For the parameters of f Orionis 
in TableQ] N* ~ 7.3 x 10 21 crrT 2 . For the line-of-sight that would apply to a neutron star emerging 
from an eclipse by f Orionis, y ~ n/2 leading to Nr ~ 1.1 x 10 22 cirT 2 , only about a factor of 
two lower the observed value in Vela X-l. The evidence, therefore, is that X-rays produced at 
R = 1 .78/?, in the wind of a supergiant companion of similar type to f Orionis are absorbed. This 
would equally be expected to apply to any intrinsic X-rays produced here and lead to the rapidly 
increasing emission measures and decreasing velocities plotted in Fig. [9] 



7. A new paradigm for the X-ray emission from O stars 

We propose that the high inertia of the plasma flow in O-star winds, caused by the weakness 
of Coulomb interactions, has important consequences for the production of X-rays in the flow. 
If shocks develop at all, they are collisionless shocks in which the fundamental role of plasma 
instabilities suggests that the local magnetic field could be the "missing variable" responsible 
for the scatter of an order-of-magnitude or more in the modest X-ray luminosities generated by 
otherwise apparently similar O stars. 

The long Coulomb collisional mean-free-path, apart from requiring a collisionless shock 
transition and slow post-shock energy exchange, also more generally limits the steepness of 
pressure gradients that can be sustained by the medium and defines a minimum size for any 
structures in the hot gas. In this case, it is unlikely that the microscopic instability in the 
line-driving mechanism will be able to steepen into the macroscopic shocks widely thought 
to generate O-star X-rays: the instability would soon be limited by the difficulty of transfer- 
ring momentum quickly enough from the unstable driven ions to the rest of the wind, caus- 
ing brea kdown of the single-fluid approx imation as discussed, for example, for thin hot-star 



winds by 



Springmann & Pauldrach 



dl992[) . Furthermore, the narrowness of the X-ray line pro- 



fil es would appear to exclude the type of positive-going velocity instabilities in the calculations 



of 



Owockietal 



(119881) . for example. 

The observed common line profile implies a common seat of X-ray emission for all the ions, 
a notion qualitatively consistent with an origin in the terminal velocity regime. While X-rays are 
probably produced throughout the wind as a routine aspect of the chaotic, supersonic flow of 
magnetized plasma, the observable volume is likely to be limited to some extent by photoelectric 
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absorption, giving rise to small systematic changes in the line profiles from ion to ion. The line 
width is mainly due to randomization of the directed kinetic energy of the flow. 

A thorough quantitative assessment of these ideas will have to await a detailed emission 
model that is expected to include, among other things, ionization, excitation and charge exchange 
averaged over the microscopic velocity distribution of ions in the post-shock gas as well as an 
energy budget with a non-thermal component. 

An inevitable consequence of the scheme proposed is the effective absence of hot electrons 
in an O-star wind, altering the physical basis of the plasma emission models concerned. In con- 
trast to either ionization by electron impact or photoionization, which together account for the 
majority of observed X-ray plasmas, single O-star spectra may be one of the clearest examples 
of a protoionized plasma. The term "protoionized" seems appropriate both for the contrast with 
photoionized and because it describes the very earliest stages of post-shock relaxation through 
which tenuous plasmas are bound to pass before electrons become hot enough to take over. 
The shock transitions themselves, though probably smaller in physical extent than the exten- 
sive shocks which span a large-scale colliding-wind flow such as WR 140, initially obey sim- 
ilar jump conditions. The distinction between spectra lies rather in the post-shock layer, in the 
amount of equilibration that takes place between ions and electrons. If the hot plasma is confined 
by magnetic fields, as in WR 140, for example, post-shock relaxation may take its long-term 
course allowing energy to be transfered from ions to electrons, which may then excite a familiar 
plasma spectrum characteristic of collisional ionization equilibrium. Otherwise, in the winds of 
single O-stars, the plasma is not in equilibrium, few or no electrons reach high temperatures and 
we observe instead the effects of the coincidence between the macroscopic terminal velocity of 
the wind and the microscopic Bohr velocity characteristic of electrons in bound atomic states. 
Interactions with the majority cool material may also be important. 

There are two simple tests of the new scheme which the data available thus far are not quite 
able to support: long, high-resolution, exposures of single O-stars should be able to establish 
two characteristics of the X-ray spectrum: first, that all the lines have much the same shape; and 
second that the continuum is weaker than that from an electron-excited plasma. 

Acknowledgements. Many thanks are due to Ton Raassen of SRON Utrecht in the Netherlands, who cast a 
critical eye over many stages of the development of the ideas discussed here. 
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Table 7. Fluxes of the X-ray lines in g Orionis using the TriLine model of Table [3] for each 
line. 



Ion 




Line ID 


Flux(10~ 5 cnr 2 s _1 ) 


Cvi 


33.734 


Lya 


40.90±1.72 




28.465 


Ly/3 


6.76+0.73 




26.990 


Lyy 


1.59±0.49 


Sxn 


36.398 




4.66+0.97 




36.573 




1.31+0.88 


S xm 


35.667 




6.96+0.94 




37.598 




3.99+1.78 




32.242 




1.93+0.68 


Nvi 


28.787 


r 


7.88+0.79 




29.084 


i 


7.48+0.85 




29.534 


f 


0.00+0.25 




24.898 


He/3 


1.35+0.73 


Nvn 


24.779 


Lyof 


17.64+0.72 




20.909 


Ly/3 


0.26+0.47 


Ovu 


21.602 


r 


96.13+2.73 




21.807 


i 


75.99+2.47 




22.101 


f 


7.69+0.96 




18.627 


He/3 


8.41+0.57 




17.768 


Hey 


1.98+0.36 




17.396 


He5 


1.60+0.36 


O vm 


18.967 


Lya 


92.43+1.73 




16.005 


Ly/3 


13.55+0.75 




15.176 


Lyy 


7.05+0.75 




14.820 


Lyd 


1.12+0.29 


Fe xvn 


15.015 




37.68+1.10 




15.262 




17.36+0.83 




15.453 




4.45+0.45 




16.778 




19.45+0.70 




17.053 




30.05+1.65 




17.098 




17.03+1.65 




12.266 




1.92+0.69 


Fe xviii 


14.208 




7.77+0.42 
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Ion 




Line ID 


FluxClO^cirrV 1 ) 




14.373 




4.03±0.36 




14.534 




2.79+0.33 




15.625 




2.91+0.39 




15.870 




1.93+0.35 




16.071 




2.96+0.61 




17.623 




2.76+0.38 


Fe xix 


13.518 




2.54+1.26 




13.795 




3.99+0.38 




15.079 




1.23+0.90 


Neix 


13.447 


r 


20.35+0.82 




13.553 


i 


11.95+1.07 




13.699 


f 


3.14+0.39 




11.547 


He/3 


2.22+0.22 




11.000 


Hey 


0.97+0.17 


Fe xx 


12.824 




0.45+0.08 




12.846 




0.45+0.08 




12.864 




0.45+0.08 


Fe xxi 


12.284 




1.27+0.29 


Nex 


12.132 


Lya 


12.02+0.45 




10.238 


Ly/3 


1.35+0.16 




9.708 


Lyy 


0.33+0.16 




9.481 


Lyd 


0.07+0.15 


Mgxi 


9.169 


r 


2.86+0.18 




9.231 


i 


1.31+0.15 




9.314 


f 


1.25+0.13 




7.851 


He/? 


0.26+0.12 


Mgxn 


8.419 


Lya 


0.55+0.13 




7.106 


Ly/3 


0.10+0.09 


Si xm 


6.648 


r 


0.98+0.12 




6.688 


i 


0.38+0.14 




6.740 


f 


0.89+0.11 


Si xiv 


6.180 


Lya 


0.27+0.11 



List of Objects 

'£ Ononis' on page| 



